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m ■ Abstract 



The scintillation efficiency (quenching factor) of fluorine and calcium nuclei 

o 



dark matter searches. The recoiling nuclei were obtained via the ^^F(n,n') and 
^'^Ca(n,n') reactions, where the neutron beam was provided via the d(d,n)^He 
reaction. The quenching factor of F and Ca nuclei in the CaF2(Eu) crystal 
^ . was found to be 11 ~ 20%, and 9 ~ 23% relative to the electron equivalent 

energy for 53 ~ 192 keV and 25 ~ 91 keV recoil energies, respectively, and 
the energy dependence was observed. The quenching factor we measured here 
is a little larger than that of previous studies, which may depend on the Eu 
doping concentration of the crystal. 
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I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs), which are major candidates for cold dark 
matter in the universe (especially in conjuction with proposed supersymmetric particles |jl|), 
interact rarely with a nucleus. When they do, the nucleus receives a recoil energy in the 10 
~ 100 keV range by elastic scattering via either the spin- independent or the spin- dependent 
interaction [^. Compared to coherent (spin-independent) scattering, experiments so far 
have far below the necessary sensitivity to exclude the dark matter candidates that have 
axial vector coupling with matter (spin coupled dark matter) |^. Recently several notable 
advantages of the spin- dependent interaction channel have been noticed @| HH [§. Calcu- 
lations by Ellis et al. suggest that fluorine provides a particularly high elastic scattering 
rate for WIMPs that have axial coupling. The rate estimation for ^^F is fairly stable against 
the calculations of nuclear matrix elements, which in some cases vary more than an order 
of magnitude for other nuclei Pj. We used fluorine-containing CaF2(Eu) scintillator. Use 
of scintillators allows fabrication of multi-kg or heavier detectors. Experiments to search 
for rare events require a large quantity of material in a low background environment. High 
statistics is especially vital for studying the annual modulation, which is a unique signature 
of halo dark matter in our galaxy. We developed a new detector system, ELEGANT VI, 
which consists of 25 CaF2(Eu) crystals totaling 7.2 kg. Details of ELEGANT VI are given 
elsewhere || 0. 

The light output of scintillators for nuclear recoil is known to be lower than that for 
electrons, particularly in the low energy region. The observed electron equivalent energy 
(pulse height) Eq is given as 

E'o = fE^R, (1.1) 

where Er is the recoil energy and / is the quenching (conversion) factor of the pulse 
height for a recoiling nucleus {€} with respect to the electron. It is known that the f- 
value depends largely on scintillators and recoiling nuclei. Although the measurements 



have been carried out by the UK (United Kingdom) Collaboration [^ |T^ and the 
BPRS(Beijing/Paris/Ronie/Saclay) Collaboration ||T3| for calcium and fluorine, these mea- 
surements are not consistent with each other. The difference may be attributed to the 
characteristics of each scintillator, such as the activator concentration, the configuration, 
and the PMT photocathode sensitivity. Therefore we studied one of the 25 crystals that are 
used for the dark matter experiment. This is essential to evaluate our detector's sensitivity 
to WIMPs. 

II. EXPERIMENTAL SETUP 

A. Pulsed-monoenergetic neutrons 

The response of the detector to the nuclear recoil of ^^F and ^°Ca was measured, where 
recoil nuclei were produced by the ^^F(n,n') and the ^°Ca(n,n') reactions, respectively. The 
experimental setup is schematically shown in Fig. |l|. A pulsed 0.75 MeV deuteron beam 
was provided by the 3.2 MV Pelletron Accelerator at the Tokyo Institute of Technology. 
The average beam current was 2 /iA at a repetition rate of 2 MHz. The ~ 3.7 MeV pulsed 
neutron beam was obtained by the d(d,n)^He reaction. The neutron flux was about 2x10^ 
/s/sr. The energy spectrum of the neutron beam was measured by a time-of-flight (TOF) 
method between the RF signal (the primary target) and a hquid scintillator placed in the 
forward direction. This indicated a spread in neutron energy of around 0.26 MeV (1 a). The 
main contribution to this spread comes from the deuteron energy loss in the primary target 
of Tid2. The Tid2 had a thickness of 0.5 mg/cm^ and was made by d2 gas electrolyzed 
on a titanium support (which also served as a deuteron beam stop). The neutron pulse 
height response of the BC501A was found to be in agreement with the calculation of V. V. 



Verbinski et al. [|14|. Details of the pulsed neutron beam are described in Ref. |T^. Neutrons 
emitted at 0° were collimated using a cylindrical paraffin collimator, 60 cm in length with a 
tapered hole along the beam axis. The hole was 1 cm in diameter near the primary target 



and 3 cm in diameter near the secondary target, the CaF2(Eu) crystaL 

B. The CaF2 detector 

The CaF2(Eu) crystal, a 45 mm cube with a europium doping content of 0.17 % by 
mass was placed at the secondary target point. It has a light guide of CaF2(pure) crystals 
coupled to photomultipliers p |]lOl- The CaF2(Eu) and CaF2(pure) crystals were made by 
BICRON. The tube is 1.5 inch Hamamatsu H3178 (quartz window). With this design the 
CaF2 detector is capable of detecting a low energy signal (0.35 keV/photoelectron). The 
whole crystal is covered with four teflon sheets (BICRON BC-642 PTFE Reflector tape) 
and an aluminized mylar ( Tokyo-met alizing 90313) to maximize the light collection. The 
linearity of the energy calibration for the CaF2(Eu) was checked using the ^'^^Am source and 
light emitting diode(LED). Typical energy resolution as measured with the ^^^Am 59.5 keV 
7 ray is about /\E / E = 25 %(FWHM). The LED (NLPB500 Nichia Chemical) is a blue 
LED and its maximum fluorescence emission- wavelength is close to that of the CaF2(Eu). 
This light is sent to both sides of the crystals through optical fibers (BCF98 BICRON). 
The energy calibration and efficiency of signals below 5 keV are difficult to test using a 
radioactive source, so we generated the pulse shape of the low energy signal with the LED 
connected to a shaper amplifier. It is verified that the signal has a similar shape to the real 
gamma ray signal by using a Flash ADC. The stability of the system was also monitored 
periodically by lighting the LED with a rate of 10 Hz during the experiment. 

C. The neutron detector and Nal detector 

The signal from the CaF2 detector was measured in coincidence with scattered neutrons 
detected by a BICRON BC501A liquid scintillator, 20 cm 0x5 cm, viewed by a Hamamatsu 
R4144 photomultiplier. This was placed 2 m from the CaF2(Eu) crystal with scattering 
angles of 30°, 40°, 50°, and 60°. Kinematics then defines the recoil energy in the target. 
Corresponding recoil energies are 53 ~ 192 keV and 25 ~ 91 keV for F and Ca, respectively. 



The liquid scintillator is sensitive to both 7-rays and neutrons. They were separated by 
using a standard pulse shape discrimination (PSD) method. The typical distribution of the 
total integrated charge over 100 ns (which depends on scintillation decay time) versus the 
initial integrated charge over 26 ns (which is proportional to the rise time) recorded by the 
ADC is shown in Fig. 0(A). The two loci correspond to the 7 and n events. The leakage 
of the 7 events to n events in the PSD selection is less than 1 % with a E^ (7 energy) > 
0.7 MeV. A heavy metal with a thickness of 10 cm was placed between the primary Tid2 
target and the neutron detector to shield neutrons and 7 rays coming directly from the Tid2 
target. Accidental coincidence due to the counting rate of the neutron detector became 
tolerable using this shielding. Accidental coincidence background was further reduced by 
the pulse height cut of 7-rays and neutrons in the neutron detector. Only events which 
produced a recoil proton pulse height in the neutron detector with a E^ (neutron energy) > 
1.2 MeV were included so as to reduce background coincidences, which come from detection 
of scattered neutrons. The neutron detector was calibrated with the '^^Y 898 and 1836 keV 
7-ray. 

Two Nal(Tl) scintillators, each with the dimensions of 85 x 95 x 95 mm^ were placed 
above and below the CaF2 scintillator to detect escaped 7-ray events. Nal(Tl) detector was 
calibrated with ^^Na 511 and 1275 keV 7-ray. The paraffin collimator prevented the Nal(Tl) 
detector from seeing the primary target position directly. 

D. Electronics 

The decay time of the CaF2 detector is quite long (~ 3/isec). The output pulses from the 
PMTs were fed into the high-gain Timing Filter Amplifier (TFA), which we developed with 
an RC circuit of integration time-constant 200 ns. The signals from the two photomultipliers 
were separately discriminated at the level well below one photoelectron. Coincidence of 
two signals from the left and right PMTs with a resolving time of 200 ns made a trigger 
signal. Pulses from two PMTs were summed and fed again into the TFA with an integration 



time-constant 1 /is. This procedure is particularly useful for a low energy signals with a 
few photons scattered during the long decay time. Then pulses above a threshold of 3 
photoelectrons occuring in coincidence with neutron events in the neutron counter were 
accepted. We developed a low noise circuit and system and were able to detect signals down 
to 3 photons under a vast amount of noise from the accelerator. 

The pulse heights from the CaF2, BC501A, and Nal(Tl) detectors and the timing infor- 
mation of the BC501A, CaF2, and the deuteron beam pulses as determined by an inductive 
pickup coil were recorded event-by-event through CAMAC-SFVME-SUN(Spark Station) on 
a disc. The trigger rate was typically about 100 cps, which contained a few cps of true 
signals. 

III. RESULTS 

Energy spectra for the CaF2 detector were obtained in coincidence with the neutrons 
selected by the PSD method. In the WIMPs search, events from the central CaF2(Eu) 



crystal are selected by setting a proper window on the roll-off ratio spectrum ||T^, which 
gives substantial reduction of backgrounds. Roll-off ratio is defined by the relative pulse 
height difference between the left and right PMT to identify the firing position in the CaF2 
detector. The light guide of CaF2(pure) crystals acts not only as a light guide but also as 
an active veto counter for the central Eu doped one. We applied the same cut used for the 
roll-off ratio. 

An energy spectrum for CaF2(Eu) at a neutron scattering angle of 60° is shown in 
Fig. 0(B). The two peaks correspond to the Ca and F recoils with energy of 91 keV and 
of 192 keV, respectively. The energy spectra were also obtained at 30°, 40°, and 50°. The 
peaks were fitted with a Gaussian to determine the mean and the width of the observed 
energy distribution. The widths of the Gaussians used in the fit include the measured energy 
resolutions of the detector, the uncertainty in the recoil energy which is caused by the spread 
in neutron beam energy, and the spread in the scattering angles. The relative intensities 
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of the two peaks versus scattering angles of neutron are consistent with calculations based 
on the known elastic scattering cross sections. It is noted here that the energy threshold of 
CaF2 detector is only a few keV, which is low enough to observe the low recoil energy. 

The fluorine peak has a tail on the high energy side in Fig. 0(B). The tail may be related 
to 7-rays from the recoiling nucleus that Compton scatter in the target scintillator depositing 
energy in addition to the recoil energy. Since the CaF2 is a low Z (atomic number) material, 
the photoelectric effect is of minor importance. ^^F has two low-lying excited levels at 109.9 
and 197.2 keV. The differential cross sections for inelastic scattering for these excited states 



are substantial according to the NNDC database [|l^. A simulation with GEANT shows 
that 21 % of 7-rays originating from these excited states interact within the crystal through 
Compton scattering. We observed that this tail was reduced by several % by requiring 
anticoincidence between the upper and lower Nal(Tl) counters. This is consistent with an 
estimate, which includes Nal's acceptance and cross sections at each angle. 

Fig. 0(A) shows the measured energy in keV for Ca and F recoils. The obtained f- 
values are also shown in Fig. ^(B) and (C) together with the previous results of the UK 



Collaboration |T2[ and the BPRS Collaboration |T^. The dominant contributions to the 
vertical error bars of the present result. Fig. ^(B) and (C), are due to the CaF2(Eu) detector 
resolution and the uncertainty in the calculated recoil energy, which originates from the 
beam energy spread due to energy loss at the primary target. 

IV. DISCUSSION 

Our result indicates that the f- value of CaF2(Eu) crystal increases with decreasing recoil 
energy, which is consistent with the result of the UK Collaboration |Tl]| [|1^. On the other 



hand, the BPRS Collaboration |T3| suggests that the f- value does not depend on the energy. 
The large f- value at low energy gives us a good sensitivity as a dark matter experiment. A 
low recoil energy is beneficial for a dark matter search. We applied two assumptions for 
the dependence of the f- value. Generally the f-value is large for particles with small dE/dX 



(electrons) |T^. We then assume the f- value depends inversely on stopping power, which is 
represented as: 

{dE/dX)0r,y\ ininiirLum— ionization 

Here {dE/dX)p(^ca) represents the stopping power of ^^F (^°Ca), for which numerical values 
are tabulated by L. C. Northcliffe and R. F. Schilling [||]. In Fig. 3(B) and (C), the f- 
values calculated by eq. (2) are compared with the present experimental data. The above 
assumption gives good agreement with the observed data. For the dark matter search we 
extrapolated the f- value to the low energy region as shown in Fig. 3(B) and (C). Typically 
dE/dX starts decreasing below the Bragg-peak, which is the energy region of interest. 
We therefore think the extrapolation to the low energy region is justified. Also shown in 
Fig. 3(B) and (C), for comparison, is a prediction from the Lindhard model \T^\. It gives a 



good agreement with ionization detectors such as Ge and Si pO| , though this is inconsistent 
with our results. 

Our f- value is a little larger than that of the other measurements. This may represent 
a dependence on the Eu activator concentration of the CaF2(Eu) crystal. We compared 
our result with the previous experimental results for CaF2(Eu) whose Eu doping content 
is known [jlT| [|l^. The observed f- values for the fluorine nucleus at a recoil energy of 100 



keV are shown in Fig. ^ There is no expression for the f-value in terms of the activator 
concentration, though a theoretical expression for the relative luminescence efficiency as a 



function of activator concentration is given E^ |23] ||23|. If we assume that we can apply it 



to the f-value also, then we obtain the curve shown in Fig. ^. Here the f-value is represented 

as, 

c(l - cY 
f « — 71 7 (4-2) 

C + (Tryl — C) 

where c is the mole fraction of the activator Eu, ar is the ratio of the capture cross section 
for the exciting energy of luminescent at a lattice trapping (non-activator) site to that at an 
activator site, and z is the number of lattice positions surrounding a given activator. The 
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parameters for the curve in Fig. 4 are 150 and 2 x 10~^ for z and a,., respectively. Our 
values for the parameters turns out to be similar to those obtained in relative luminescence 
efficiency measurement of Nal crystal, which are 60 and lxlO~^ for z and 0"^, respectively 
P3|| . However, there are 3 parameters including normalization for 3 data points, so more 
data is necessary. This assumption indicates that the f-value of our crystal would have an 
activator concentration close to that giving the maximum f-value. Since there is no clear 
foundation for using eq. (3) for the f-value, further studies are needed. 

V. CONCLUSION 

In summary, we observed recoils of calcium and fluorine nuclei in the range 25 ~ 91 keV 
and 53 ~ 192 keV in a CaF2(Eu) crystal, respectively. The scintillation efficiency was found 
to be 11 ~ 20%, and 9 ~ 23% for F and Ca, respectively. These data serve to calibrate a 
CaF2(Eu) dark matter detector whose results are discussed in Ref. [ p^ . 
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FIGURES 
FIG. 1. Schematic top- view of neutron scattering apparatus. The neutron collimator was made 

from a cylindrical paraffin, 60 cm in length with a tapered hole along the beam axis. The hole was 

1 cm in diameter near the primary Tid2 target and 3 cm in diameter near the secondary target, 

the CaF2(Eu) crystal. Each scintillator was viewed by a photomultiplier tube (PMT). 

FIG. 2. (A): The typical distribution of the total integrated charge over 100 ns (which depends 
on scintillation decay time) versus the initial integrated charge over 26 ns (which is proportional to 
the rise time) recorded by the ADC. (B): An energy spectrum of CaF2(Eu) at a neutron scattering 
angle of 60°. The two peaks correspond to Ca and F recoils with energy of 91 keV and 192 keV, 
respectively. 

FIG. 3. (A): The measured energy in keV for Ca and F recoils in CaF2(Eu) calibrated relative 
to electons (7-rays) at each energy, versus the corresponding calculated recoil energy. The errors 
are statistical only. 

(B) and (C): The f- values (quenching factors) versus recoil energy in CaF2(Eu) for (B) Ca recoils 
and (C) F recoils. The crosses are the data from the present experiment, diamonds and stars are 



data points from the UK Collaboration |12| and the BPRS Collaboration |13|, respectively. The 



f- values are not plotted by the BPRS Collaboration [13|, so we converted by taking only their 
vertical (statistical) error bar into account. The systematic errors are also included for the UK 
Collaboration and ours. The dashed line is the extrapolation of the f-values calculated by eq. (2) 



and Ref. ||18[ (squares). The solid curve is the prediction from the Lindhard model [19|. 



FIG. 4. The f-values as a function of Eu concentration (mole fraction) for fluorine nucleus in 
CaF2(Eu). Experimental data from Ref. [jll|(open circles) and this work(full circle) are shown, 
and the errors are only statistical. The dashed line is obtained from eq. (3). 
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